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Abstract 

Einstein-Infeld-Hoffmann method is used to solve the problem of motion 
of two bodies when the equations of general relativity are of the generalized 
form: they have been reduced to a form invariant under conformal trans- 
formations. It is proved that not only metric degrees of freedom, but also 
derivatives of vector A a appearing in the generalized equations can exert in- 
fluence on the motion of bodies in a certain space-time domain. This influence 
can account for the recently observed anomalous acceleration of spacecrafts 
Pioneer 10, Pioneer 11. The impact of vector A a on the motion of bodies is 
interpreted as a consequence of viscosity in geometrodynamic continuum. 
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1. Introduction 



The analysis of the paths of spacecrafts Pioneer 10, Pioneer 11 in refs. [1], 
[2] suggests that in range 

R = (30 4- 100) AU = (0.45 4 1.5) • 10 15 cm (1) 
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the spacecrafts are subject to an anomalous acceleration component, which 
has come to be denoted by ap in the literature. In magnitude the acceleration 
ap is 

a P = - (8.74 ± 1.33) ■ 1(T 8 cm/s 2 . (2) 

The minus sign indicates that the direction of the acceleration ap is close to 
the direction towards the Sun. 

From the viewpoint of the following discussion of the issue of the accel- 
eration ap it seems reasonable to make two estimations. First, we should 
compare (0) to acceleration a^ gained by spacecrafts from the Sun in the 
Newtonian approximation. Having taken the mass of the Sun equal to 
M — 2 ■ 10 33 g, gravitational constant G = 6.67 ■ 10~ 8 cm 3 /g ■ s 2 , we ob- 
tain for distance R = 50 AU = 7.5 ■ 10 14 cm 

GM (6.67 -10- 8 cm 3 / g ■ s 2 ) (2 ■ 10 33 g) 9 „ 7in _ 4 /2 ^ 
aN = l¥ = 5.625 ■ 1(P = 2 - 37 ' 10 ° m/S - (3) 

Second, for the system (the Sun + spacecraft) we should determine small- 
ness parameter A, which is on the order of magnitude of the ratio between 
the characteristic velocity of the craft relative to the Sun and light velocity. 
For the above value of R: 



IGM 



(6.67 • 10" 8 cm 3 /g ■ s 2 ) ■ (2 • 10 33 g) . _,, 



c 2 R \ (3 • 10 10 cm/s) 2 ■ 1.5 • 10 15 cm 



10" 5 . (4) 



The above estimations show that the ap is not predictable by the general 
relativity either in the Newtonian or post-Newtonian (PN) approximation. 
If fact, the Einstein-Infeld-Hoffmann (EIH) corrections to acceleration differ 
from accelerations by a value of the order of ~ A 2 , while relation {ap/a^l 
is close to A 1 . The corrections to the acceleration in the PN approximation 
coincide with the distance and the ap is independent of distance. 

It is our view that from the standpoint of explanation of the nature of ap 
the fact is essential that product He, where H is Hubble constant, is equal 
to 

He = 6.9 • 10~ 8 cm/s 2 . (5) 

which on the order of magnitude is close to ap (it is accepted that H = 
2.3- 10~ 18 1/s), that is 

ap ~ —He. (6) 
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The closeness between ap and He is noted both in ref. [1] and in many other 
papers. In this connection many people put forward the hypothesis that the 
additional acceleration ap owes its origin to the cosmological expansion of 
the Universe. However, no acceptable theoretical construction implementing 
this hypothesis has been proposed. 

This paper makes an attempt to consider the motion of bodies in systems 
like (the Sun + spacecraft) within the framework of a so-called conformal 
geometrodynamics (CGD). CGD is a theoretical scheme based on equations 
that are a minimum conformally invariant generalization of Einstein equa- 
tions for empty space. The CGD equations are derived in [3] (without A 
term) and in [4] (with A term). They are: 

R a p — ~9apR = —2A a Ap — g a pA 2 — 2g a pA u v + A a .p + Ap. a + A • g a p. (7) 

Here A a , A are the vector and scalar fields that are attributes of the Rieman- 
nian variety dynamics. The form of equations (|7j) is preserved in transfor- 
mations 

g a /3 -> g a p ■ e 2<J , A a > A a a. a , A -> A • e~ 2<T , (8) 

where a is an arbitrary scalar function of coordinates. These transforma- 
tions pertain to the category of those, which effect the conformal mapping 
of Riemannian spaces and are termed conformal transformations (see, e.g., 
[5])- 

To find solutions to equations (J7J), coordinate conditions and gauge must 
be specified. For the coordinate conditions the well-known de Donder condi- 
tions will be used generalized so that they acquire the form invariant under 
transformations (JSJ). The conformally invariant form of the de Donder con- 
ditions is: 

-T^(V=gg M ) +24,-^(7=^ = 0. (9) 

V-9 ' a V-9 

For the gauge condition this paper uses condition 

A = A = Const. (10) 

When (jl(Jj) is valid, vector A a automatically satisfies Lorentz condition: 

A* = 0. (11) 

As will be demonstrated in this paper, the consideration of the body 
motion dynamics in a system like (the Sun + spacecraft) allows realization 
in a seminal manner of the hypothesis that the additional acceleration ap 
follows naturally from the CGD equations. 



3 



2. Equations and conditions in the EIH scheme used 



We will solve the problem of finding the equations of translational motion 
of two point bodies in the lower orders of approximation in parameter Xv/c, 
where v is the characteristic velocity of relative motion of the bodies. We 
will be solving with the Einstein-Infeld-Hoffmann (EIH) method in the form 
presented in refs. [8], [9]. When comparing specific expressions of this paper 
with the relevant expressions in refs. [8], [9] it should be kept in mind that 
this paper uses signature ( — h ++) and for the coordinate condition the de 
Donder condition in form @; the other notations are the same as those used 
in refs. [8], [9]. 

Assume that in our approximations A exerts no influence at all on the 
body motion and can be excluded from what follows 1 . 

In the EIH method, the arrangement of orders of smallness of j a p, A a , 
where 

la/3 = h a p - -r) a /3 (lf V hp V ) , h a /3 = Qa.fi ~ Vafi, 

Va/3 = diag [—1, 1, 1, 1] is the metric tensor of plane background space, de- 
pends on parameters of the bodies (masses, velocities) and distances between 
them. As applied to the system (the Sun + spacecraft) at distances ((H), this 
consideration adopts the following order of smallness arrangement: 



7oo = 
7ofc = 

7mn 

A = 
A k = 



7oo + 7oo + 7 o + -, 

2 4 5 

:7ofc + 7 fc + 7 fc + ..., 

3 4 5 

= 7 mn + 7 mn + 7 mn + 
4 5 6 

A + Ao + 

3 4 

Ak + A k + ••• • 

4 5 



(12) 



As it will follow from the following, the above order of smallness arrangement 
rules out the influence of field A a on spacecraft motion in the Newtonian 
approximation. 

Equations (J7J) lead to the following equations for j a p (only those equations 
are written out which are used in what follows): 



00: A z 



0. 



(13) 



1 For this reason hereafter there will be no danger to confuse function A (x)with small- 
ness parameter A. 
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Ok; A ; 



-- (Al ok ) + - f-7oo,o + 7 0s ,s 

2 V 3 / 2 V 2 ' 1 3' 



0. 



(14) 



[mn; A 4 
i 



> i I 7 Om, n 7 ms.sn I 
1 ' 3 1 4 



■7 



On, m + 7 ns,sm 
3 1 4 



+ I 7 00 7 00,mn~ 



"mn ( 7 00, — 2 7 0s, s + 7 

\2 1 1 3 1 4 , y . 

+ S 7 00,m 7 00,n — l^ran 7 00 ( A 7 00 ) — Tannin 7 00,s 7 
^ 8 2 2 8 2 V 2 / iD 2 2 

J 2 V 4 



A 7 mn 
4 



(15) 



00,; 



• 7 oo, o + 7 os, 



= 0. 

A o,fc- 



00; A z 



L 1 3 3 /. 

-- A 7 oo + 77 7 pqm + — 7 00 , s 7 00)S + - 7 oo A 7 oo 

9/i 2 4 y y 16 2 2 8 2 V 2 



[mn; A 5 ] 
1 



•> ) 1 ' 7 0m, n + 7 ms,sn + — 7 On, m + 7 
z 1 \ 4 ' 1 5 / V 4 1 5 



$mn ( 2 T 0s, s + 7 
V 4 l 5 

A m,n A n,m 2(5 m 



pq.pq 



(16) 

= 0. 
(17) 

(18) 



^4 o,o + A i,i 

4 5 



Coordinate conditions Q are of form (in square brackets: C.C. =Coordinate 
Condition): 

aC.;0;A 31 



C.C; k; A 4 j 

GC.;0;A 4 



7om = 7 00 ,o+2Ao- 

3 2 l 3 



7 fci / — 7 ofc, o — t 7 oo 7 oo,fc- 

4 3 1 4 2 2 ' 



7om = 2 ^o- 

4 4 



(19) 

(20) 
(21) 



C.C; 0; A 



1 1 1 

7oz,z — 7 oo, o — -7op7oo,p — t7oo7oo,o + -7oo7os,s + 2v4o- 

5 4 1 23 2 42 2 22 3 5 



C.C.;A;;A £ 



7ofc,o-7w,i + 2A* = 0. 

4 l 5 5 



(22) 
(23) 



Among all conditions (jllj) . only the gauge conditions of the lower approx- 
imations are sufficient (in square brackets: G.C =Gauge Condition): 



G.C; A 



A i,i 

4 



0. 



(24) 



G.C.-X 



A 0,0 + A 1,1 = 0. 

4 5 



(25) 



From equations (0) it follows that functions AoAo, Ak, Ak should sat- 

3 4 4 5 



isfy the following equations: 



0; A 3 




-Aio 

3 


= 


0; A 4 " 




-AAo 

4 


= 


'k- A 4 " 




4 


= 


~k; A 5 " 




-AA fc 

5 


= 



(26) 

(27) 
(28) 
(29) 



Thus, functions Ao, Aq, A hi Ak should be harmonic functions. 

3 4 4 5 



3. Centrally symmetric static solution to CGD equa- 
tions 

The centrally symmetric static (CSS) solution to the CGD equations is pre- 
sented in different forms in several papers (see [6], [7]). Here the CSS solution 
will be presented in a form convenient for comparison to the Schwarzschild 
solution. 

In the CSS problem the squared interval is 



ds 2 



-exp (7) ■ dt 2 + exp (a) ■ dz 2 + exp (f3) ■ (d0 2 + sin 2 ^ ■ dyj 2 ) . 



In this case the gauge vector can have as few as two nonzero, radial coordinate 
dependent components: 

A a = (0,/,O,O). 

Having made use of the gauge transformation with radial variable dependent 
function a = o~(z), set A\ — f equal to zero. Then replace the radial 
coordinate so that condition goo^n = — 1 be met. This transformation does 
not result in appearance of ^-component of field A a . Upon the above two 
transformations, without loss of generality, the metric can be written in the 
form 



ds 2 



-exp (7) • dt 2 + exp (-7) • dz 2 + exp (/3) • (d9 2 + sin 2 6> • dip 2 
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and vector A a in the form 



0,0,0). 



Four equations are obtained for four functions 7, (3, 0, A. 
Equation G° = T ° 



exp (7) • (3" + f {(3'f + \i(3' - exp (-/?) = 3exp (-7) • 2 + A. 



Equation = Tl 



exp (7) • [I (/?') 2 + IV/3'J " exp (-/?) = exp (-7) • 2 + A. 
Equation G 2 . = T| 

exp (7) ■ [±/3" + I + |V + I ( 7 f + I^Vl = exp (-7) ■ + A. 



Equation Gq = Tq 



= 0'- 7 >. 



It turns out that there are three solution types. We are concerned about 
that type, which contains the Schwarzschild solution as a special case. Here 
is the solution type with the procedure itself of its finding being omitted. 

= p ■ exp (7) 
exp (J3) = A ■ sh 2 (p z + o ) 

+ B ■ [p z ■ cth (p z + a ) - 1] + b p ■ cth (p z + a ) 



ex P (7) = ^ 
X(z)=B p 2 

(30) 

The zero-subscripted quantities are integration constants. 

With appropriately chosen constants a solution of type (J30|) can approx- 
imate the Schwarzschild solution as closely as is wished in a certain range of 
the radial variable. To obtain the approximation, set the following in ()3()j1 : 



B = 0, a = 0, A pl = 1. 

Then for e 7 : 

e 7 = 1 + 6 p • cth (p z) . 
Having introduced notation b = — r , = 1/A write the CSS solution 



1 ^ 4-V 

1 • cth 

L 



LJ ' 
7 



^ = L 2 ■ sh 2 



(31) 
(32) 

(33) 



L [ L \L 
In the range of the radial variable 



(r /L) < (z/L) « 1 



(35) 



expressions (I33|) . (J34)) assume the following form: 




(36) 




z 



(37) 



Expressions (J3fij) are the same as the associated expressions in the Schwarz- 
schild solution. In so doing ro has the meaning of gravitational radius. As 
for (|37[). there is no analog of this expression in the Schwarzschild solution. 
In expressions ()36|). ()37|). two facts seem essential. First, in the range of 
radial variable ()35|). for its description the CGD equation solution requires 
not only dimensional constant ro, but also one more dimensional constant L. 
Second, in the above-mentioned range of radial variable the principal term 
in the expansion of function is constant 1/L. 

4. Finding functions y a p, A a with the EIH method 

For us it is convenient first to find the equations of translational motion of two 
particles under the assumption that the particle masses are commensurable 
in magnitude, that is, to obtain results similar to those obtained in refs. [8], 
[9], however not for the equations of general relativity, but for equations (J7J). 
Then from the found equations of motion the equations of motion can be 
obtained straightforwardly for the system, in which one of the particles is a 
trial particle. 

The general approach to the solution of our problem is that components 
7 Q/ 3 will be represented in form j a p = 7a/3 + Sj a p. Here 7<*/3 are expressions 
that are obtained in the framework of general relativity for two point particles 
of masses M,m without inclusion of additional terms of A a . On the other 
hand, £70/3 is an addition determined entirely by the A a . 

The structure of all quantities j a /3, A a is essentially predetermined in the 
scheme under discussion by expressions for 700 and A in the lower orders 
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of smallness. In our case the expressions for 7 00 , 7 00 should coincide with 

2 3 

expressions 

Mm . . 

7oo = 4— + 4-, 7oo = 0, 38 
2 n r 2 3 

that is with those expressions, which in the case of general relativity lead to 
the well-known equations of translational motion in the PN approximation. 

As for v4o, certain considerations regarding its choice follow from approxi- 
mate expressions (j3T?j) . (|3T|) for the exact CSS solution to the CGD equations. 
Clear that the expansion of A should begin with the constant 1/L, where 
L is the dimensional parameter used to determine asymptotic A at long 
distances. When considering multi-particle problems, the expansion of the 
exact one-particle solution does not allow detection of the presence of parti- 
cle velocity containing terms. Appearance of this type of terms cannot also 
be excluded in the expression for Aq, when a system of at least two particles 
is considered. Note that a similar situation arises, for example, when one 
attempts to find the Kerr solution with the EIH method [9] or tries to deter- 
mine radiation friction force in the Dirac-Lorentz equation with this method 
[10]. 

Reasoning from the above considerations, write the expression for A , 
which could be composed of the quantities present in the problem. The 
desired expression: 

- should begin with constant 1/L, 

- can contain particle velocities r]k in the next order of smallness, 

- should be, according to PBj) . (J2ZJ), a harmonic function, 

- if it contains a part depending on coordinates, then the part should be 
decreasing with distance from the system, 

- should be symmetric about replacement of parameters of one particle 
by those of the other. 

It turns out that all of these requirements are satisfied, if Aq is represented 

as 

A ° = i + i 0= L +K - (m + M)L\— + 4 } (39) 

and R/L is assumed to be of the third order of smallness. Here k is a 
constant numerical coefficient, fi = mM/ (m + M) is scaled mass of a system 
of two particles. In this consideration multiplier fi/ (m + M) is assumed to 
be of the order of unity. The need for the introduction of k to the A 
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construction is due to the lack of uniqueness that is inherent in the EIH 
method in construction of approximations to the exact solution with using it 
(for more details, see [9], [10]). The considerations on selection of the value 
of the constant k will be presented later. 

Ak is found from gauge conditions From equation and 

condition (J24j) it follows that A k — 0. Write the expression for A^ in the 

4 

vicinity of the first particle alone; in so doing make use of the equation of 
motion of particle in the Newtonian approximation. 



A k = A k + A k 

4 5 



mR k n {RiRittk FPlX&fa 



(m+M)L | Rr 1 r\ rf 



(40) 



From here on the term with two arrows (1^2) denotes the expression de- 
rived from the written one through replacement of parameters of one particle 
by those of the other. 

The substitution of expression for A Q into equations (JHJ), (fTtjj) with 
taking into account (|19jl. (J2"T|) yields the following form of the expression for 

Sj ok = Slo k + Slo k = ^- k .^ + (1^2). (41) 

3 4 oL [m + M) L n 

Found expressions (p^). (f4*Tj) are sufficient to determine the form of 
the correction to the Newtonian equation of particle motion. These correc- 
tions are obtained by integration of equation (fTH|) over infinitesimal-radius 
spheres surrounding the point particles. Thus we obtain the correction to 
the equations of motion for the first particle. 

Write equation (fTHjl with having rearranged all the terms to the left side. 

ms,sn + -7 On, m + 7 ns.sm 
2 I V 4 1 5 ) V 4 1 5 ' 

A T rran ^mn ( 2 7 Os, s 7 pg,pg ) f (^2) 
5 \ 4 1 5 / J 



-<4 m,n A n,m 20j-, 
5 5 



A 0.0 + A M 

4 ' 5 



0. 



In the integration of this equation the total contribution from the terms 
containing second derivatives of components £7 mn is zero, as they form curl 
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combinations. The contribution from the following combination is zero for 
the same reason: 

~ \ 8 mn 7 Os, s "Ton.Om 
2 IV 4 1 4 1 



The contribution from 



+25 m „ (-Ao,o 



is zero by virtue of the gauge condition. As a result, to find the correction 
to the equation of motion, the contribution from the following terms must 
be calculated: 

— r 7 Om, n + -5 mn 7 Os, s — A m,n n.m- 
2 4 l 2 4 l 5 5 

This correction will therewith enter into the equation of motion 

mM 

2M£ k = -jp-R k + 3 fe (43) 

as addition 

^ k = \n i \r\ l ofc ' 5" + 2^ fcn l ° s ' i s _ 4 fc,n ~ 4 n ' k ] dSn ' ^ 

The contributions of separate terms to surface integral f)44j) are: 



4-7T 



f ||4n7os, o s | rfs n = ^ § ^5 kn Ao,o\ ds n 

A k ,n- An,k^ds n = K- (-§ • f ■ ^ - | ■ ^ • (RiRi)i k 

Substitution of (ft5j) -(|47 ]l into (jSJ) gives: 

Sfc = K . . 1^ + 2k ■ L— ■ (Rtft) 4. (48) 

(m + M) L R (m + M) L \ 1 l )^ k ^ > 



(46) 



(47) 
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As a result, from (|43|) . (J44|) . ()48|1 we obtain the following correction to 
acceleration S^k- 



^ k = K - 2( m TM)M -lk + K - ( m + M)ML -( RlRl )^ (49) 



In terms of derivatives with respect to regular time and masses in regular 
units of measurement equation (J4TJJ) is written as: 

d 2 5£ k _ fJ.m c 2 R k fxc 2 ( dRA d^k , , 

dt 2 ~ K ' 2 (m + M)M' LR (m + M) MLG ' \ 1 dt J dt ' 1 j 

Assume that constant L is of cosmological origin and is related with 
Hubble constant H as 

L = c/H. (51) 
Having substituted (fBHj) into (jBTljl . we obtain: 



d 5£k Kfxm TT Rk ( dRi \ <i£, 



dt 2 2 (m + M) M C R (m + M) MG ' ° ' { 1 dt J dt ' 1 ' 

Expression (jB^j) is just the correction to the Newtonian expression for 
acceleration which follows from the CGD equations under the assumptions 
specified at the beginning of this section. 

5. Equations of motion of spacecrafts Pioneer 10, 11 

Write equations (|H2~|) in the special case, where mass of the first particle is 
much less than that of the second particle, that is where the first particle can 
be considered clS £L trial particle. In this case 

^ = *.cH R ±-^.cH.( Ri §l)^. (53) 
dt 2 2 R mG \ dt j dt K ' 

It is interesting to note that the first term in the right-hand side of (I53j) 
depends neither on masses of the bodies nor on the distance between them 
whatsoever. 

From (|53J) it follows that in the trial particle motion along circumference 
the correction from the second term in the right-hand side of (|53|) is zero, the 
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correction to the Newtonian expression for the acceleration is determined by 
the first term. 

d 2 8£ k k R k 

~W = 2 ' CH ll- (M) 
Consider motion of the trial particle along radius. We will omit subscript 
"1" in writing the radial vector values. In this case we obtain for the radial 
component of acceleration from (|53jl : 

d ^ = --cH-^.cH.R( d 4 d f). (55) 
dt 2 2 mG \dtdtj V ; 

Write the squared radial velocity with making use of the law of conser- 
vation of energy, that is as 

2E 2Gm 
+ 




M R 



We arrive at: 



dt 2 2 mG V M R J 2 mMG v ; 

Consider the spacecrafts as trial particles. Total energy E of the space- 
craft leaving the system with minimum kinetic energy is close to zero. As- 
suming E = in (57), we obtain 

^ = -«.?.cff. (57) 
dt 2 2 K ' 

To make formula (joT|) consistent with formula (jUJ), it should be assumed 
that the coefficient k is close to unity. More precisely, 

if K=l then ^f = -cH 1 

if k=1, then ™f = -l. c H J " [ ' 

We arrive at the conclusion that by fitting one constant multiplier CGD 
allows the observed anomalous component of the acceleration of spaceships 
Pioneer 10, 11 to be described. The available experimental data on the 
spacecrafts Pioneer 10, 11 can be described by formula (|57|) with k ranging 
from k = 2/3 to K = 1. 
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6. Conclusion 



From formulas (J57|) . (|58|) it follows that as the spacecraft moves away from 
the Sun along radius: 

1. The additional acceleration is directed oppositely to the spacecraft 
direction and is constant in magnitude. 

2. As long as the spacecraft can be viewed as a trial body, the acceleration 
is independent of the spacecraft characteristics and is of universal nature. 

It should be emphasized that all the consideration is valid only when the 
initially made assumptions are fulfilled. Below they are mentioned in the 
explicit form. 

First , the space-time dynamics is described by the generalized equations 
of general relativity - conformal geometrodynamics equations (J7|). 

Second , the principal term in the expansion of Ao is a constant in the 
space-time domain under consideration and the next term is a quantity of the 
fourth order of smallness. These assumptions are nontrivial, as in this manner 
cosmological-origin quantity L = c/H is introduced to the construction of 
the Aq to determine asymptotic expressions for the components of vector A a 
that satisfy equations (|IJ). In the general relativity there is no vector A a , so 
L cannot be introduced naturally in it. 

In the space-time domains, where the above two assumptions are not 
fulfilled, first, expression (}39|) forAo can have a higher order of smallness, 
second, terms of another type (for example, pole terms) can become principal 
in the construction of Aq. 

Expression (p)3|) for the additional acceleration of the trial body admits a 
direct experimental verification. Thus, it predicts: 

- universal nature of the additional acceleration for all bodies moving 
along circular orbit [only the first term in the right-hand side of (J53|) remains] , 

- possibility to change the acceleration direction in the motion along ra- 
dius (in moving away - towards the Sun, in approaching - away from the 
Sun). 

In conclusion note that in equations (jo^j) . (jo^j) the additional acceleration 
owes its appearance either to the terms including derivatives of A a vector 
components or to terms that vanish in disappearance of this derivative type. 
The thermodynamic analysis performed, in particular, in ref. [7] suggests 
that these terms determine the geometrodynamic continuum viscous stress 
tensor. Therefore we conclude that the additional acceleration owes its origin 
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to the geometrodynamic continuum viscosity. 

The appearance of the anomalous component (provided that it is not 
explained by non-gravitational effects) raises a number of basic issues. For 
example, as applied to the general relativity, this is the issue of energy con- 
servation law feasibility degree and cosmological expansion effect on the ex- 
periments within the Solar system. 

Spacecraft experiments can be undertaken in the future in order to arrive 
at answers to fundamental questions of the space-time theory. It is not 
improbable that the results of this paper will prove useful in development of 
programs of investigations in these experiments. 

This work was carried out under financial support by the International 
Science and Technology Center (ISTC Project #1655). 
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